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We present a quenched calculation of meson correlators with the overlap fermion 
at very small masses 2.6-13 MeV. In this region the pion Compton wavelength is 
larger than our lattice size (L ~ 1.23 fm) and the system is in the so-called e-regime 
of chiral perturbation theory. We found that the scalar and pseudo-scalar correlators 
are precisely approximated by a few hundred low-lying fermion eigenmodes in this 
regime, whereas axial-vector correlator receives significant contributions from higher 
eigenmodes. We also measure the disconnected pseudo-scalar correlator, which is well 
saturated with the low-lying modes. Matching these lattice data with the one-loop 
expressions for the correlators in quenched chiral perturbation theory, we evaluate the 
decay constant F„ and the chiral condensate E as well as the parameters wig and a, 
which describe the artifacts of the quenched approximation. 
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§1. Introduction 



Chiral perturbation theory (ChPT) provides a systematic method to calculate low energy 
dynamics of QCD, though it contains unknown parameters at each order of the expansion 
in pion mass squared m 2 and momentum squared p\. At the lowest order these parameters 
are the pion decay constant F w and the chiral condensate E, and there are 10 other low 
energy constants in the next-to- leading order. 1 ) It is one of the important tasks of lattice 
QCD to calculate these low energy constants non-perturbatively from the first principles. In 
the standard approach, however, this is very demanding because one has to work on large 
enough lattices satisfying m^L 3> 1 to avoid possible finite size effects. Taking the chiral 
limit — > and the continuum limit while satisfying this condition will be prohibitively 
time consuming even with today's fastest supercomputers, especially when the quarks are 
treated dynamically. 

In the so-called e-regime, 2 )~ 4 ) where the linear extent of the space-time box is smaller than 
the pion Compton wave length L <C l/m w (but larger than the QCD scale 1/vIqcd <C L, 
which assures that the pion can be treated as a point particle and other heavier hadrons are 
decoupled. ), the chiral Lagrangian is still applicable except that the expansion parameter 
is given by e 2 ~ m n / 'A ~ p^/A 2 , where A is a cutoff scale of the chiral Lagrangian roughly 
around 1 GeV. An important observation is that the low energy constants in the chiral 
Lagrangian are defined at the cutoff scale commonly for both the standard and e-regimes. 
Therefore, one can determine the low energy constants in the e-regime and use them in the 
standard ChPT. In this way one can avoid the problem of the chiral limit while keeping the 
large lattice volume. Analytic calculation of meson correlation functions in the e-regime is 
available for both quenched and unquenched theories. 5 )~ 7 ) 

In the study of the chiral regime of lattice QCD the chiral symmetry plays an essential 
role. First of all, one must treat pions near the massless limit, which appear as a result of the 
spontaneous chiral symmetry breaking. Furthermore, it is known that the effects of fermion 
zero modes become important in the e-regime 8 ) and the correlators depend strongly on the 
topological charge of the background gauge field. 5 )' 6 ) The lattice fermion formulations pre- 
serving chiral symmetry 9 ) (by satisfying the Ginsparg- Wilson relation 10 )) is now commonly 
used (but only in the quenched approximation). In this work we use the Neuberger's overlap- 
Dirac operator. 11 )' 12 ) With this formulation there is no fundamental problem to approach 
the chiral limit as required in the study of the e-regime, but the computational cost to invert 
the overlap-Dirac operator increases for small quark masses. 

In the chiral regime the meson correlators are largely affected by the low-lying fermion 
modes especially by the chiral zero modes. In this work we explicitly study the effects of 
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such low-lying modes using the eigenmode decomposition of the fermion propagator. In 
the quenched approximation we find that the connected scalar and pseudo-scalar meson 
correlators are reproduced to 98-99.9% accuracy (depending on channels) with only 200 
lowest-lying eigenmodes on a 10 3 x 20 at small quark masses (2.6 MeV < m < 13 MeV). 
Such saturation was previously found in Refs. 13), 14), but it should be better in the e- 
regime (In our study, L satisfies m n L ~ 0.6 and AqcdL ~ 3.). An advantage of such 
eigenmode decomposition is that the meson correlators can be averaged over space-time 
points without much extra computational costs. The statistical fluctuations originating from 
the local bumps of zero-mode wave function are suppressed by the space-time averaging and 
thus we can avoid the large noise as found in Ref. 15). The low-mode averaging was also 
used in recent studies. 16 ^' 17 - ) 

Matching our numerical data of the axial-vector, scalar and pseudo-scalar correlators 
with the quenched ChPT (QChPT) expressions, we extract the leading order low-energy 
constants £ and F w as well as the parameters appearing due to the quenched artifact, a and 
m . The axial- vector correlator is most sensitive to F n while £ is precisely determined by the 
connected scalar and pseudo-scalar correlators. We also investigate the chiral condensates 
and the disconnected (hairpin) correlators for the pseudo-scalar channel. In general we find 
good agreement between the lattice data and the QChPT predictions for topological charge 
\Q\ =0 and 1 sectors, but the larger topological sectors deviate significantly, which may 
suggest breakdown of the e-expansion in the QChPT at large \Q\. 

The axial- vector correlator has already been calculated in the recent works. 15 )' 17 ) Bi- 
etenholz et al. 15 ' worked at a relatively larger quark mass (~ 21 MeV) and found that the 
correlator is well fitted with the QChPT formula for a large enough lattice size (L > 1.1 fm). 
They also pointed out that the signal is very noisy at \Q\ =0. Using the low- mode averaging 
technique, 16 )' 17 ) Giusti et al. pushed the quark mass down to 10 MeV and found an encour- 
aging agreement of F n measured in the e- regime as that from the standard measurement. 
Their result is F n = 102(4) MeV in the e- regime. There is another interesting work by the 
same authors, 18 - ) who investigated the divergent (as ~ 1/m 2 ) contributions of zero mode in 
the massless limit and matched them with the theoretical expectations from QChPT. 

This paper is organized as follows. In Section we review the quenched chiral per- 
turbation theory (QChPT) in the e-regime following the discussions of Damgaard et al. 5 ^ 
We describe the details of our simulation in Section |3] and study the low-lying eigenmode 
dominance in Section 0] In Section |S] we present our results for the chiral condensate and 
meson correlators and the comparison with the QChPT. Conclusions are given in Section El 
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§2. Quenched chiral perturbation theory in the e-regime 

In this section, we briefly review quenched chiral perturbation theory (QChPT) in the e- 
regime 5 ** and summarize the relevant formulae for our analysis of meson correlation functions. 
The partition function of QChPT with N v valence quarks is written as 

Z(9,M) = J dUexp y- J d*x£ M (x)J , (2-1) 

where the Lagrangian C e M is given by 

£mW = IfStT&Uixy^Uix)) - ^fstr(U e U(x) + U{x)- l U^) 

2 

777 C\ 

+^$(x) 2 + —dMx)dMx) (2-2) 

at the leading order of m 2 and p 2 expansion. The field variable U(x) is integrated over a sub- 
manifold of the super-group Gl(N v \N v ), the maximally symmetric Riemannian sub-manifold, 
which is characterized by a matrix of form 



U = ( „ , Ae U(N V ), D e Gl(N v )/U{N v ), (2-3) 




and Grassmannian N v x N v matrices B and C. Str denotes the super-trace. The mass term 
corresponds to the choice of mass matrix M = (ml v +ml v ) with I v and I v the identity matrix 
in the fermion-fermion and boson-boson blocks respectively. The effect of CP violating 9 term 
enters through Ug = exp(i9/N v )I Nv + I Nv . In the quenched approximation the singlet field 
${x) = ^Str[— iln U(x)] does not decouple, and the couplings Wq and a are introduced. 19 -* 
The e-regime 2 -' -4 ) is realized when the quark mass is small enough that the pion Compton 
wavelength ~ l/m^ is larger than the linear extent of the space-time L. The systematic 
expansion is then reorganized and the expansion parameter is given by e 2 ~ m 7T /A7cF n ~ 
l/(LF n ) 2 . Unlike the standard ChPT the zero mode of pion gives important contribution 
and one must explicitly integrate out the constant mode of U(x). This is done by writing as 

U(x)=U ex pi y^- (2-4) 

and integrating over the constant mode Uq. One obtains the partition function at a fixed 
topological charge Q by Fourier transforming (j2-l|) . 



Zq(M) = -!- [ +W d9e WQ Z(9,M) 
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'"~ V Str(U + U 



+ j d'x (-istr(0 M £0M0 - ^(StrO 2 - ^(^StrO 2 ) + 0(e 



2 

1 (2-5) 



where dUo denotes the Haar measure of the maximally Riemannian sub-manifold of Gl(N v \N v ). 
The topological charge distributes as Gaussian with variance 

V ~ 2N C [Z ' b) 

in the quenched theory, which is an exact equation for V — > oo. It shows a good contrast 
with the full theory, for which (Q 2 ) = mEV/Nf is expected for Nf flavors. 

We note that in the quenched approximation the Gaussian approximation of the Fourier 
transform in 6 is justified only for small topological charge that satisfies \Q\ <C (Q 2 )- (See 
Appendix for the details.) Therefore, all the results shown below are valid only for small 
\Q\. We investigate how this breakdown of the effective theory occurs using the lattice data. 

In the following we consider N v = 1 and 2 as we are interested in the system with two 
light quarks. All the results are obtained by the perturbation of £ fields and the exact 
integration over zero mode Uq, which can be written in terms of the Bessel functions. 

At the tree-level the scalar condensate is given as 

- (#> Q = EqQm) = Eh(i m (h)k m m + i M+1 (n)K\ Qh M) + (2-7) 

with n = mEV . I\q\{^) and K\q\(h) denote the modified Bessel functions. The \i dependence 
of Eq(/j) is shown in Figured] Near the massless limit it asymptotically behaves as E\q\ (//) — > 
E\Q\/fi for \Q\ > 0. One- loop correction does not change its functional form 20 ) 

£q~ 1oOP ()«) = S e«fi'(I\Q\{^')K m (ii') + I m+1 (jj,')K\Q h i(ji')) + r eff i^i = E Q (n'), (2-8) 
but the parameters \i and E are shifted to // and E e g: 

fi' = mE cS V, (2-9) 
^ SE ( 1+ n#H0)+*X0) j ( , 10) 

Here, parameters (5(0) and A(0) are ultraviolet divergent tadpole integrals, 

1 e ipx 

p^O y 

1 e ipx 
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Fig. 1. Eq{h) in different topological sectors at = 270MeV. 

which need to be renormalized. In our analysis they are to be determined by matching with 
lattice data. 

Let us define the flavor-singlet meson operators 

S°(x)=^(x)I Nv ^(x), (2-13) 
P\x)=^(x)i lb I Nv ^x). (2-14) 

Adding these operators to the QCD Lagrangian as source terms 

C ^ £ + s(x)S°(x) +p(x)P°(x), (2-15) 

corresponds to a substitution 

M - M + s(x)I Nv + ip{x)I Nv (2-16) 

in the effective theory. The two-point correlation functions of them are obtained by dif- 
ferentiating the generating functional with respect to s(x) and p(x). To 0(e 2 ) the results 
are 



(S°(x)S°(0)) Q = C» + 



(P°(x)P°(0)) Q = CI 



E 2 
2F 2 

E 2 
2F 2 



{m 2 Q G(x) + aA(x)) - A(x 



a + + a_ 



^KG(x) + aA(x)) - A(x) a+ + a 2 +A 



where 



a+ = 4 
a_ = 4 



jJL E 



V 2 



A* 2 J 



(2-17) 
(2-18) 

(2-19) 
(2-20) 
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and the constant terms are written as 



C° P 



y eff 1-loop _ ^2 

— a + - A, eS 



y 2 

^eff 1-loop _ v-,2 



E, 



eff 



Q 2 

jl' 1 



i W) q 2 

H' E cS /i' 2 



Note that the prime denotes the derivative with respect to ji, 



\Q\ 



For flavor non-singlet mesons, we need a N v — 2 super-group integral, which 
described by the Bessel functions. The non-singlet operators are given by 

S a (x)=i,(x)(r a /2)I Nv ij(x), 
P a {x)=^{x){T a /2)i l5 I Nv ^{x). 

with the Pauli matrices r a . To 0(e 2 ) the two-point functions are given by 



>23) 
is also 



(S a (x)S a (0)) Q = C a s + 
(P a (x)P a (0)} Q = C P - 



E 2 
2F 2 



E 



2 r 



2F 2 



-r(m 2 G(x) + aA{x)) - A{x)b_ 
^±(m 2 G{x) + aA(x)) - A(x)b+ 



where 



2 1 + 



o: 

V 2 



fi 2 



c + = 2 
c_ = -2 



E 

/i £ 



and 



£ra = ^eff ( ^Q(^')\ 



2 



E, 



eff 



2 \ /iTeflf 

For the flavor non-singlet axial-vector current 

A»(x)=^(x)(r°/2)i 7/l 7 6 ^(x) 



(2-34) 
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the correlator is obtained as 6 -* 



(A a o (x)A a o (0)) Q = -y- 2mZ Q ^)A{x) 



(2-35) 



Here, an important observation is that the axial-current correlator does not involve the 
parameters related to the quenched artifact, i.e. tuq and a. We also note that the constant 
term is proportional to rather than E. Therefore, this channel is suitable for an extraction 
of F n , whereas for the pseudo-scalar and scalar correlators F w appears only in a coefficient 
of A[x) and G(x) terms. 

In the lattice calculation we measure the correlators with zero spatial momentum pro- 
jection. It is therefore convenient to define 

h 1 (\t/T\) = ^jd 3 xA(x)= 1 - 



\t\ 
T 



1 

12 



(2-36) 



h2(\t/T\) = ~ [d*xG(x) = ± 



J^2 



M-iV-i. 

T ) 30 



(2-37) 



In the e-regime the correlators do not behave as the usual exponential fall-off exp(— Mt) with 
the mass gap M to be expected in the large volume. Instead, it becomes a simple quadratic 
function for the single pole A(x) and quartic function for the double pole G(x) integral. 

As these expressions show, the meson correlators in the e-regime are quite sensitive to 
the topological charge and the fermion mass. Hence they provide a good testing ground for 
lattice simulations in the e-regime. Furthermore the parameters F w , E, m and a can be 
extracted from the fitting of these correlators. The parameter £ always appears associated 
with the quark mass m. It makes sense because only the combination mU is renormalization 
scale and scheme independent. The numbers we extract for S in the following analysis should 
be understood as a result in the lattice regularization at a scale 1/a. To relate them with 
the conventional scheme such as the MS scheme requires perturbative or non-perturbative 
matching, which is beyond the scope of this paper. 

§3. Lattice simulations 



We generate the gauge link variables at (3 — 5.85 in the quenched approximation on a 
10 3 x 20 lattice. The lattice spacing is a = 0.123 fm, which is obtained from the Sommer 
scale r = 0.5 fm using an interpolation formula given in Ref. 21). The linear extent of the 
lattice is then about 1.23 fm. We employ the overlap-Dirac operator defined by 



D rn = (l-—)D + m, 
D = -(l + -y 5 8ga(H w )), 



(3d) 
(3-2) 
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Table I. Number of configurations in each topological sector. 

\Q\ 12 3 

# of confs. 20 45 44 24 



with the kernel H\y built with the Wilson-Dirac operator D\y, 

H w = 75(aD w - 1 - s). (3-3) 

The parameter s controls the negative mass given to aDyy and we choose s = 0.6 at (3 — 5.85 
to minimize the number of low- lying mode in H\y. The symbol sgn denotes a sign function 
of the large sparse Hermite matrix H\y, and a is defined as a = a/(l + s). This overlap-Dirac 
operator satisfies the Ginsparg- Wilson relation 

l5 D + D l5 = aD l5 D (3-4) 

exactly, and the 75-hermiticity = 75-D75 is also satisfied. In the practical implementation 
we approximate the sign function sign(ifjy) using the Chebyshev polynomial of degree 100- 
200 after subtracting 60 lowest-lying eigenmodes of sign(H\y) exactly. The error of the sign 
function is then 10~ 12 level and safely neglected in our numerical results. 

One of the essential points of our work is to use the eigenmode decomposition of the 
fermion propagator. For this purpose we calculate 100 lowest (but non-zero) eigenvalues of 
P±DP± and their eigenfunctions as well as zero-mode eigenfunctions of P T DP T for negative 
(and positive) topological charge. We use the numerical package ARPACK, 22 ) which imple- 
ments the implicit restarted Arnoldi method. The chiral projection operator P± = (l±7 5 )/2 
is applied in order to reduce the rank of the matrix. The eigenvalues and eigenvectors of the 
original matrix D can be reconstructed from those for the chirally projected operators. We 
thus obtain 200+ \Q\ eigenmodes of D for each gauge configuration. Note that these 200+ \Q\ 
eigenvalues cover more than 15% of the circle in the complex space of the eigenvalues of D 
as Figure El shows. The topological charge is obtained from the number of zero-modes and 
their chirality. The number of configurations for each topological sector is given in Table HI 
We analyze the gauge configurations of \Q\ < 3 sectors. 

When the exact inverse of the overlap-Dirac operator is needed, we use the techniques 
described in Ref. 23). For a given source vector 77, we solve the equation 

D m ip = r] (3-5) 

by separating the left and right handed components as ip = P-tp + P+ip and solving two 
equations 

P^ = (P„DiD m P„y l P_Di V , (3-6) 
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0.5 1 1.5 2 

Fig. 2. Lowest 202 eigenvalues of the overlap-Dirac operator at beta=5.85 on a 10 3 20 lattice of 
topological charge Q = —2. The eigenvalues cover a ir/3 arc of the circle. 

P + ij = {P+D m P + )-\P + r) - P +J D m P_V), (3-7) 

consecutively. (The above equations apply for positive Q and the same procedure applies 
with a replacement P + <-> P_ for negative Q.) We use the conjugate gradient (CG) algorithm 
to invert the chirally projected matrices with the low-mode preconditioning with 20 lowest 
eigenmodes. With this low-mode preconditioning we gain about one order of magnitude 
speed-up of the CG solver (we need only 20-40 iterations for each P±ip.) for our smallest 
quark mass 0.0016, which corresponds to 2.6 MeV in the physical unit. 

§4. Low-mode dominance for the meson correlators 

The inverse of the overlap-Dirac operator D can be decomposed into the contributions 
from each eigenmode with eigenvalue Aj and eigenvector Vi(x) as 

-^low 

D-\x,y) = n - /on T < x >l{v) + AD m \x,y). (4-1) 
z — ' (1 — am i) A* + m 

i=l 

Here, the eigenmode decomposition is done incompletely, and the sum is truncated at some 
cutoff iV"i ow , which we set iV low = 200 + \Q\. The additional term AD^{x,y) represents the 
contribution from higher eigenmodes. 

We expect that the low energy physics is dominated by the low-lying eigenmodes. Near 
the massless limit the lowest-lying eigenmodes play a dominant role as they are enhanced by 
a 1/ Xi factor, and especially the zero modes give divergent contribution in the quenched ap- 
proximation. (In the unquenched case their occurrence is suppressed by the fermion determi- 
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nant.) Sensitivity to the gauge field topology mostly comes from these low-lying eigenmodes. 
Therefore, the eigenmode decomposition (|4TJ1 is expected to give a good approximation for 
low energy physics even if we ignore the higher mode contribution AD^-(x, y). The cutoff 
iVi ow must be large enough for such approximation to cover the all relevant low-lying modes, 
and it depends on the pion mass and physical volume of the system. For massive fermions 
the eigenmodes below Aj < m become equally important, and we need larger iVi ow than in 
the massless limit. 

On the other hand, the short distance physics could be affected by the higher eigenmodes. 
It can be seen by looking at the defining equation D m (x, z)D^(z, y) = Xa=T Vi{x)v\(y) + 
D m (x, z)AD^ n 1 (z,y) = S(x,y) at x = y. Since the sum J2f=T v i( x ) v l( x ) approaches mono- 
tonically to one when iVi ow is sent to the size of the matrix -/V max , the contribution from the 
remaining term AD^{x,y) becomes significant for N\ ow much smaller than iV max . However, 
such a short distance correlation described by AD^ (x, y) should be insensitive to the gauge 
field topology. 

4.1. Connected correlators 

First, let us consider the "connected" meson correlators 

Mi(t)Mi(o) conn = ^^r^(f,t)vSr^(o,o) COIin 

X 

= - E tr { r * D m@ > *5 0, O^D-^O, 0; x , t)) , (4-2) 

X 

where Mj denotes local operator corresponding to pseudo-scalar (PS), scalar (S), vector (V) 
and axial- vector (AV) currents, and denotes the corresponding gamma matrix, i~p S = 
r s = 1, r v = ij Q , and r AV = «757o- 

We calculate these correlators in two ways: one is an exact calculation with the conjugate 
gradient (CG) method and the other is the low-mode approximation 

200+IQI 

where the higher eigenmode contributions are neglected. Comparison of them at m = 0.0016, 
0.0048, 0.008 (< 13 MeV) is shown in Table ITTl where maximal numerical difference between 
the two correlators in the region 7 < t < 13 is listed for each operator and topological 
sector. For the scalar and pseudo-scalar mesons, the approximation (J4-3)) does work well to 
98-99.9% accuracy, while the lowest 200 + \Q\ modes are not enough to reproduce the vector 
and axial vector correlators. Figure H3 shows the pseudo-scalar and axial- vector correlators 
at m = 0.008 and \Q\ = 1. We observe very good agreement for the pseudo-scalar correlator 
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Table II. Comparison of the low-mode approximated "connected" correlators with their exact 
ones. The approximation uses the quark propagator (|4-3|) with 200 + \Q\ low- lying modes at 
m = 0.0016, 0.0048, 0.008. The maximum deviation in 7 < \t\ < 13 is shown. The number of 
configurations is given in Table EJ 



m = 0.008 (~ 13 MeV) 


101 =o 


101 = 1 


101 =2 


101 =3 


scalar 


2.01% 


1.49% 


0.46% 


0.44% 


pseudo-scalar 


1.19% 


0.48% 


0.28% 


0.22% 


vector 


41.6% 


258 % 


259% 


94.6% 


1 j 

axial-vector 


28.4% 


a r\ t- fry 

40.5% 


24.3% 


24.2% 


m = 0.0048 (~ 7.7 MeV) 


101 =o 


\Q\ = 1 


\Q\ =2 


101 =3 


scalar 


1.57% 


0.65% 


0.16% 


0.18% 


pseudo-scalar 


1.06% 


0.33% 


0.10% 


0.08% 


vector 


28.9% 


198 % 


110% 


104% 


axial-vector 


26.2% 


38.6% 


19.7% 


18.3% 


m = 0.0016 (~ 2.6 MeV) 


101 =0 


101 = 1 


101=2 


101=3 


scalar 


1.41% 


0.08% 


0.03% 


0.03% 


pseudo-scalar 


0.95% 


0.06% 


0.04% 


0.02% 


vector 


22.5% 


149 % 


157% 


166% 


axial-vector 


21.5% 


32.6% 


20.4% 


34.7% 



for wide region of t (from 3 to 17), but the axial- vector case is much worse. This is probably 
because the axial- vector correlator is smaller in magnitude by a factor of 0(m), and therefore 
small fluctuation of eigenvectors is enhanced. 24 -* Figure 0] shows how the low-lying mode 
contributions saturate to the full correlators. We find that about 100 lowest eigenmodes 
suffice to approximate the full correlator to a very good accuracy. The plot is shown for 
m = 0.008, but the saturation becomes even better for smaller quark masses. 

An advantage of the low-mode approximation (|4-3J) is that D^(x,y) at any x and y is 
obtained without performing the CG inversion, so that one can easily average the source 
point over the space-time 



^^(f,t)^(0,0) ) ^_^/^^(f,* + *o)^(f ,*o)) • (4-4) 



This so-called low-mode averaging dramatically reduces the fluctuation of the low-lying 
modes as shown in Figure El which is also reported in Ref. 17). In practice we average 
only over {L/2) 3 x (T/2) lattice points where the site index is an even number for each 
direction. Note that even after the low-mode averaging, the error from the truncation of the 
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Fig. 3. Pseudo-scalar (left) and axial-vector (right) correlators at m = 0.008 and \Q\ = 1. 
low-mode-approximated correlator (|4 ■!-{[) is compared with the corresponding exact one. 
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Fig. 4. Saturation of the connected pseudo-scalar correlator for one sample configuration at m = 
0.008. Two plots are shown for \Q\ = (left) and \Q\ = 1 (right). The rightmost points around 
120,000 correspond to the exact correlator obtained by the CG method. 



higher modes is negligible compared to the statistical error ~ 15% in Q = sector and ~ 
5% in Q ^ sectors in the range 7 < t < 13. 

4.2. Chiral condensates 

Consider the low-mode contribution to the scalar and pseudo-scalar condensates for a 
fixed topological charge Q 

(^(x)f = -(trD- 1 (x,x)f 

'Mow 1 \ \ Q 

E n - /oux v i( x >i(x) + AD m \x,x) ) , (4-5) 
^— j' (1 — am/2)\i + m J 

(^(x)) Q = -(tr 75J D m 1 (x,a;)) Q 
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Fig. 5. The pseudo-scalar correlator at m = 0.008 and \Q\ = 1. Filled symbols denote the data 
with the low-mode averaging, while open symbols are not averaged. 



= - ( tr (g (1 - am/2)X, + m -wWWM + T^'fe «)) ) ■ <«) 

For the scalar condensates, it is known that AD^(x,x) includes unwanted additive ultravi- 
olet divergences 25 -* 

(trD m \x, x))Q = —^-- 3 + C 2 ™ + Cl— + S Q {n') + 0{e% (4-7) 
(1 + s)a A a 1 a 

where and C\ are unknown constants. The first term comes from the modification of the 
chiral symmetry in the Ginsparg- Wilson relation, 10 ) 7 5 D _1 (a;,?/) + D~ x {x,y)^ = a8 x ^ y) at 
x = y. If iVi ow is large enough, the higher mode contribution AD^(x, x) should be insensitive 
to the link variables U^(y) separated large enough from x and thus to the global structure 
of the gauge field configuration, such as the topological charge. We, therefore, expect that 
such contribution vanishes in the difference between different topological sectors, 

-((^(x)) Q - (4-8) 

In other words this difference should be well described only by the low-lying eigenmodes. In 
fact we observe such a low- mode saturation as shown in Figures E] and while the individual 
condensate — is n °t saturated. 

For the pseudo-scalar condensate, we do not have to take care of the higher modes 
because the condensate is determined only by the zero modes and the contributions from 
other eigenmodes cancel because of the orthogonality among different eigenvectors. As shown 
in Figure |Hl our data with iV low = 200 + \Q\ low-modes perfectly agree with the theoretical 
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4.3. Disconnected correlators 



Here we discuss the expectation value of the "disconnected" diagrams for a fixed topo- 
logical charge. As in the case of the connected diagram, we expect 

disc 

£ tr i r i D m& *; £ *)) tr ( r ^m(0, 0; 0, 0) 
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Fig. 8. Pseudo-scalar condensates approximated with iV low = 200 + \Q\ low-modes. The lines 
represent the expectation (|4-9|) . 
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^ V 7 r 
tr (i-am/: 
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/2)Ai + m 



«*(0)«J(0) 



+ 2<tr(r 4 Z\D- 1 (x,x))>' 



x ( tr 
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^ E (1 - am/: 



/2)Ai + m 



«i(0K(0) 




+ 



(tr (/vaitffo *)))'" 



(4-10) 



where a; = (af, i). We assume that higher modes' contribution does not have correlation with 
any local operator 0(y) separated enough from x, i.e. 



(AD-\x,x)0(y)f (AD-\x,x))' x (O(y)Y 



(4-11) 



where the expectation value (■••)' represents insensitivity to the topological charge. We 
also use the translational invariance (0(x)) = (O(0)). Unlike in the "connected" case, we 
cannot check the low-mode dominance by explicitly computing the exact correlators because 
the numerical cost is too expensive. However, for the pseudo-scalar disconnected diagram 
(tr7 5 Z\D~ 1 (x, x))' vanishes because AD^^x, x) does not contain the zero modes. 

In fact, as Figure El shows, we find good saturation with the lowest 200 eigenmodes for 
the pseudo-scalar disconnected correlators. Similar results were also obtained previously in 
the study of the rj' propagator with the Wilson fermion 26 ) and with the overlap fermion. 27 ) 
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§5. Extraction of the low energy constants from the meson correlators 

5.1. F n from the axial-vector correlator 

The axial-vector current correlator (|2-H5jl is most sensitive to F n and not contaminated 
by the parameters mo and a. The problem is, however, that the low-mode dominance does 
not hold, and we have to solve the quark propagators exactly. Hence, the statistical noise 
could become a problem as we cannot average the source point over space-time. Our strategy 
is to treat data at different quark masses and different topologies at the same time to reduce 
such large statistical errors. 

We use the local axial-current A a ^(x) = ip(x)'~f5'~f fl (T a /2)ip{x) , constructed from the over- 
lap fermion field ip(x). Since it is not the conserved current corresponding to the lattice 
chiral symmetry, (finite) renormalization is needed to relate it to the continuum axial-vector 
current. To calculate the Za factor non-perturbatively, we follow the method applied in 
Refs. 28), 29). Namely, we calculate 

p( ) " £^(z,t)P"(o,o)) ' (54) 

where Vo denotes a symmetric lattice derivative. The pseudo-scalar density P a (x) must be 
the chirally improved one ^(r a /2)7 5 (l — ^D)ip associated with the exact chiral symmetry 
on the lattice. For the on-shell matrix elements such as the one considered here, one can use 
the equation of motion to replace the aD term by —am /{I — am/2), which is negligible for 
our quark masses. We therefore use the local operator for P a (x). 

The ratio (|5-lj) turns out to be insensitive to topology as shown in Figure Fitting 
the average of R p (t) over all topological sectors with a constant in the range 7 < t < 13, we 
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Fig. 11. ap as a function of the bare quark mass. Dashed line is a result of the quadratic fit: 
ap = 0.00001 + 1.390(ma) - 0.19(ma) 2 . 



obtain the results for 

_ a(V^(x)P»(0)) 2ma 

ap{ma) = (p*(x)p°(o)) " + 0(a } (5 ' 2) 

at four quark masses m = 0.0016, 0.0048, 0.008, 0.0160, which are shown in Figure ITT1 With 
a quadratic fit we obtain 

ap(ma) = 0.00001(2) + 1.390(14)(ma) - 0.19(74)(ma) 2 . (5-3) 

The constant term is perfectly consistent with zero and we can extract Za from the linear 
term as Za = 1.439(15) , which is consistent with the value Za = 1.448(4) reported in 
Ref. 30) which was done with the same (3 and s. 

We now compare the renormalized axial-vector correlation function with the QChPT 
result 

2Z 2 X>o(£,t)Ao(0,0)> ( 2 = 2 + 2m27|q|0*)r/i 1 (|t/r|)^ . (5-4) 
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From the constant term we can determine F n , while E has to be extracted from the small 
t dependence. In Ref. 15) it is reported that the correlators suffer from large statistical 
fluctuation at \Q\ = and other topological sectors are insensitive to E. We observed 
similarly large statistical noise as shown in Figure IT21 but it turned out that two-parameter 
(F n and E) fitting does work well when we treat the data of different topology and fermion 
masses simultaneously 

As Figure fT2"l shows, our data at m — 0.0016, 0.0048, 0.008 in the \Q\ < 1 sectors are 
well described by the QChPT formula (|5-4|) . A simultaneous fit in the range 7 < t < 13 
yields F w = 98.3(8.3) MeV and E 1 ^ = 259(50) MeV with x 2 /dof = 0.19. These results 
are stable under a change of the fit range: (F w , E 1/3 ) = (98.8(8.3) MeV, 261(47) MeV) and 
(99.2(8.3) MeV, 261(45) MeV) for 8 < t < 12 and 9 < t < 11, respectively. The result for F n 
is in agreement with that of the previous work, 17 ) 102(4) MeV. The authors of 15) quoted a 
slightly larger value ~ 130 MeV. 

The correlators at \Q\ = 2 do not agree well with the above fit parameters as shown in 
Figure Uni As discussed in Section El it may indicate that the topological sector \Q\ = 2 is 
already too large to apply the QChPT in the e-regime. 

5.2. E, E c g and a from connected S and PS correlators 

As discussed in the previous section, the scalar and pseudo-scalar connected correlators 
are approximated rather precisely only with the lowest 200+ \Q\ eigen-modes at small quark 
masses (m = 2.6-13 MeV). In this way we measure the scalar (pseudo-scalar) correlators 

(S(t)f = -2 £(1 + s) 2 (S%x)S%0))? ow _ modes , (5-5) 
(V{t))* ^2^(1 + s) 2 (P 3 (x)P%0))? ow _ modcs , (5-6) 

for the topological sectors 1 < \Q\ < 3 (0 < \Q\ < 3) at m = 0.0016, 0.0032, 0.0048, 0.0064, 
and 0.008, in which the error from higher mode truncation is estimated to be only < 1 %, 
which can be ignored compared to statistical errors. We take an average of the source point 
over {L/2f x (T/2) lattice sites. 

To fit our data with the one- loop QChPT formulae ()2-26|) and (j2-27|) . we have to determine 
five parameters in these formulae: F n , E, E E g, uIq and a. Since the sensitivity for F v is 
weak with these correlators, we use the jackknife samples of F n obtained through the axial- 
vector current correlator, F n = 98.3(8.3) MeV. Unfortunately, there still remain too many 
parameters to fit with QChPT expressions. Therefore, to determine wig we use the relation 
()2-6|) and input the value of topological susceptibility x = (Q 2 )/V from a recent work 31 ) as 
r oX — 0.059(3). It gives m = 940(80)(23) MeV, where the second error reflects the error of 



19 




0.003 



0.0025 



0.002 



0.0015 



0.001 



0.0005 



qChPT F=98MeV,£=(259MeV) J 



i' 



10 

t 



15 



20 



O 

CO 



O 

CO 

o 



< 

V 



0.003 



0.0025 



0.002 



0.0015 



0.001 



0.0005 



0.003 



0.0025 



0.002 



0.0015 



0.001 



0.0005 



qChPT F=98MeV,E=(259MeV) J 



^^••l-i-i-i-i-i-i-ff-rri" 1 ' 



10 

t 



15 



qChPT F=98MeV,E=(259MeV) J 



10 

t 



15 



20 



20 



Fig. 12. Axial-vector current correlators at m = 0.0016 (left) and 0.008 (right) for \Q\ 



and 1 (bottom). The dashed lines are the result of simultaneous fitting of the data for Q 
and 1 at m = 0.0016, 0.0048, 0.008 in the region 7 < t < 13. 
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13. Axial- vector current correlators at m = 0.0016 (left) and 0.008 (right) for \Q\ =2. The 
dashed lines are the QChPT prediction with parameters determined through \Q\ =0 and 1 
sectors. 



20 



O 



(f) 

V 



lattice Q=1 i — i- 
latticeQ=2 :— -*- 

qChPTQ=1 - 

qChPTQ=2 




10 12 14 16 18 20 
t 



o 



lattice Q=0 
lattice Q=1 
lattice Q=2 
qChPT Q=0 
qChPT Q=1 
qChPT Q=2 



s S3 s m so 



10 12 14 16 18 20 



t 



to 
O 



V) 

v 



E 

as 

o 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 



-0.2 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 



-0.2 



lattice Q=1 i — i — 
lattice Q=2 v—x— 

qChPTQ=1 - 

qChPTQ=2 




* w, m m * m 




10 12 14 16 18 20 



t 



lattice Q=0 
lattice Q=1 
lattice Q=2 
qChPT Q=0 
qChPT Q=1 
qChPT Q=2 



H H g H B o. 



10 12 14 16 18 20 



t 



Fig. 



14. Scalar (top) and pseudo-scalar (bottom) correlators at m = 0.0032 (left) and 0.0064 
(right). The dotted lines are the fit results with all available mass parameters m = 0.0016, 
0.0032, 0.0064 and 0.008 for topological charges and 1. 



roX- 

Using these inputs, we fit the correlators (|5-5j) and (|5-6|) in the range 7 < t < 13 at 
different Q and m simultaneously. Figure El shows the correlators with fit curves. For 
\Q\ < 1, the data at all available quark masses m = 0.0016, 0.0032, 0.0064 and 0.008 are 
fitted well, and we actually obtain x 2 /dof ~ 0.7. (Note that the correlations between different 
i's, m's and channels (PS and S) are not taken into account.) Our fit results are S 1 ^ 3 = 
257 ± 14 ± 00 MeV, which is consistent with Ref. 32), Z^ 3 = 271 ± 12 ± 00 MeV, and a 
= —4.5 ± 1.2 ± 0.2, where the first error is the statistical error and the second one is from 
uncertainty of (Q 2 )- They are insensitive to the choice of the fit range. The central values 
vary only slightly (for instance, ± 1 MeV for E 1 ^ 3 and S^ 3 ) by choosing shorter fitting 
ranges 8 < t < 12 and 9 < t < 11. 

From the ratio E e g/U we can identify the size of the NLO correction in the e expansion. 
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To the one-loop level, it is written as 

= 1 + j^r W<5(0) + aA(0)) , (5-7) 

where the parameters (5(0) and 4(0) are regularization dependent. For this ratio we obtain 
1.163(59), which indicates that the e expansion is actually converging. 

We obtain large negative value for a, which is also reported in Refs. 33) and 34). These 
results contradict with a previous precise calculation, 35 ) which obtained a small value a 
= 0.03(3). If we instead assume a = 0, and fit as a free parameter, we obtain F n = 
136.9(5.3) MeV and £ and E e s are almost unchanged. (Detailed numbers are summarized 
in Table IIHI ) Therefore, there is an apparent inconsistency in the determination of F v 
between the axial- vector and (pseudo-) scalar correlators if a ~ is assumed. A possible 
cause is that \Q\ = 1 is not small enough to derive the partition function Eq. (j2-5|) (See 
Appendix). Eq. ()2-6|) may also have a systematic error due to finite V as well as finite a. 

The data at higher topological charge, \Q\ = 2, are also plotted in Figure fT4l They do 
not quite agree with expectations from the QChPT shown by dashed curves in the plots. A 
simultaneous fit with all the data including \Q\ =0, 1 and 2 gives a bad x 2 /dof (~ 12). This 
problem of higher topological charge also happens for the axial- vector correlator as discussed 
in the previous subsection. 

5.3. Scalar condensate 

The free parameter in the scalar condensate is E e g as seen in (|2-8|) and (|2-9|) . To avoid the 
problem of ultraviolet divergence we compare a difference — ((ipip)® — (ipip) °) with the QChPT 
result Eq(h') — Eq=q(h'). We use the low-mode approximation with 200+ \Q\ eigenmodes, 
and the low-mode averaging is done as in the meson correlators. Figure El shows the 
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difference as a function of quark mass for \Q\ = 1, 2 and 3. We find that the lattice data 
agree remarkably well with the QChPT expectation with E e ff = 271(12) MeV as determined 
from the (pseudo-) scalar connected correlator even in higher topological sectors. In fact, if 
we fit the scalar condensate with E e ^ as a free parameter, we obtain 256(14) MeV which is 
consistent with the result above. 

5.4. Disconnected PS correlators 

We also measure the disconnected pseudo-scalar correlators, which is made possible with 
the low-mode approximation. The saturation with 200+ \Q\ lowest-lying mode is quite good 
for the pseudo-scalar channel as discussed in Section 14.31 

In the QChPT it is written as 



{V d (t)) Q = J d 3 x (1 + s) 2 (2P 3 (x)P 3 (0) - P°(x)P°(0)) Q 



d A x (I + s 



,2 



C d P - ^ + aA(x)) - e + A(x 



(5-8) 



where 



C p = Su*> ( 5 - 9 ) 
d+ = -4 (l + ^ , (5-10) 

e+ =-2(Y^Y-^#y (5-n) 



E J iiE 

In Figure El lattice data for topological sectors \Q\ — 0-3 are shown at two representative 
quark masses m = 0.0032 and 0.0064. The QChPT predictions are plotted with the parame- 
ters determined through the axial- vector and (pseudo-) scalar connected correlators: S 1 ^ 3 = 
257 MeV, F n = 98.3 MeV, m = 940 MeV, and a = -4.5. We observe that the agreement 
is marginal, though the correlator's magnitude and shape are qualitatively well described. 
Instead, if we fit the disconnected correlator with E and a as free parameters while fixing 
F w and m to the same value, we obtain E 1 ^ 3 = 227(32) MeV and a = —3.5(1.2), which 
are statistically consistent with those input numbers. Therefore, we conclude that both 
the connected and disconnected correlators are consistently described by the QChPT in the 
e-regime. Details of the fit results are listed in Table ITTT1 



§6. Conclusions 

In the e-regime of ChPT, the meson correlators are largely affected by the fermion zero- 
mode, and thus by the topological charge of background gauge field. This expectation from 
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0.0064. The curves represent the results of QChPT with E 1 ^ = 257 MeV, F n = 98.3 MeV, 
m = 940 MeV and a = -4.5. 



the effective field theory is explicitly confirmed by a first-principles calculation of lattice QCD 
using the overlap Dirac operator, with which we can preserve the exact chiral symmetry at 
finite lattice spacing. 

To reach the e-regime we need small quark masses to satisfy the condition M^L < 1. This 
can be achieved by using the eigenmode decomposition of the fermion propagator. Then, 
the connected scalar and pseudo-scalar correlators are precisely reproduced by using only 
200 low-lying eigenmodes. This number would be unchanged even when we decrease the 
lattice spacing, as far as the physical volume is kept fixed to ~ (1.2 fm) 3 x (2.4 fm), since 
the small eigenvalue distribution depends only on a combination \EV in the e-regime. The 
a — > limit affects higher modes only, which are irrelevant to the low energy dynamics. 
For those connected meson correlators the small quark mass regimes are reached without 
extra computational costs, and we can also employ the low-mode averaging technique to 
substantially reduce the statistical noise due to near-zero mode contributions. For the axial- 
vector current correlator, on the other hand, the saturation by low-lying modes is much 
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Table III. Summary of the fitting results. The first column denotes the topological sectors used 
in the fit. The values in [■ ■ ■ ] are input parameters. The first error is statistical. The second 
and third errors reflect the uncertainty in the input parameters, (Q 2 ) and F n , respectively. 



correlators 


F w (MeV) 


27V3 (MeV) 
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worse and we had to treat them exactly using the (costly) CG solver. We also investigate 
the disconnected pseudo-scalar correlator using the low-lying mode approximation. The 
disconnected diagrams are usually very expensive as they need many fermion inversions, 
but with this approximation they are obtained without extra costs. We confirm that the 
disconnected pseudo-scalar correlator is well saturated by 200 low-lying eigenmodes for our 
lattice. 

Remarkable Q and /j, dependences of the quenched ChPT are well reproduced by our 
lattice calculation at f3 — 5.85 on a 10 3 x20 lattice. Then we are able to extract some of the 
low energy constants: F w , E, and a. The last two describe the artifact of the quenched 
approximation. Fitting our data for meson correlators simultaneously at different quark 
masses and topological charges, we obtain F w = 98.3(8.3) MeV, E 1/3 = 257(14) (00) MeV 
(E 1 ^ 3 = 271(12)(00) MeV), m = 940(80)(23) MeV, and a = -4.5(1.2)(0.2), from the con- 
nected correlators with \Q\ < 1. In these numerical results the second error reflects the error 
of {Q 2 ). We also obtain consistent results from disconnected pseudo-scalar correlator and 
the chiral condensate. 

Despite such remarkable success of QChPT, we also find problems. First, the correlators 
in \Q\ > 2 sectors are not well fitted by the parameters determined at \Q\ < 1. This 
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may indicate a problem of the e-expansion in QChPT, which comes from the fact that the 
partition function at a fixed topology can be justified only for \Q\ <C (Q 2 )- Since (Q 2 ) is 
proportional to V, this condition is relaxed on larger lattices. It is therefore interesting to 
extend our study to larger volumes. Secondly, the numerical result for a is rather large and 
negative. Such a large value may pose a question on the validity of the partition function 
()2-5|) . Previous lattice calculations, e.g., 35 ** indicated the value consistent with zero and our 
result is clearly inconsistent with them. If we assume a = 0, then our data prefer larger 
value of F n ~ 130 MeV, which contradicts with the result of the axial- vector correlator. Our 
calculations are, of course, not free from other systematic errors due to higher order terms 
in QChPT, finite lattice spacing, etc. However, because of the exact chiral symmetry the 
finite lattice spacing error does not spoil the consistency with (Q)ChPT, while the extracted 
parameters are contaminated. 

Once these problems are (positively) solved, the lattice simulation in the e-regime could 
become a strong alternative to the conventional large volume (or large quark mass) simula- 
tions. A clear advantage is that the small enough quark masses can be reached and there is 
practically no question on the applicability of ChPT. Obviously, the dynamical simulations 
in the e-regime are most desirable. With the overlap fermion we do not see any fundamen- 
tal problem, since the overlap-Dirac operator is well-conditioned even in the massless limit, 
though it is numerically too expensive on the current-generation machines. 
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Appendix 



In this appendix we derive the partition function at a fixed topological charge (|2-5jl . It 
is obtained from the partition function (j2-l|) by Fourier transforming in 9. 



Z Q (M) = — I d9e t0Q Z(e,M) 



2tt 



d6 / dUod£($detU' ) Q exp 
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where we use 



U = U e^ F « 
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U' Q = UgU 
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^Str(-<lnUS)=<Po + ^. 



In the last line of (jA-ljl . we perform 6 1 integral as a Gaussian; 
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(A-6) 



(A-7) 



where (Q 2 ) = Vm^F 2 /2N C and 9' = 9 — \/2^o/ i 7 ^- To justify this Gaussian integral, we need 
a condition \Q\/(Q 2 ) <C 1, otherwise the integral (|A-7J) should depend on $q, which means 
that $o and 9 can not be treated independently and the partition function Eq. ([2-5j) is not 
valid. On our lattice with (Q 2 ) = 4.34(22), this condition becomes questionable for \Q\ ^ 0. 
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